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Abstract 

The small-XBj limit of deep inelastic scattering is related to the high-energy limit of the 
forward Compton amplitude in a familiar way. We show that the analytic continuation 
of this amplitude in the energy variable is calculable from a matrix element in Euclidean 
field theory. This matrix element can be written as a Euclidean functional integral in 
an effective field theory. Its effective Lagrangian has a simple expression in terms of 
the original Lagrangian. The functional integral expression obtained can, at least in 
principle, be evaluated using genuinely non-perturbative methods, e.g., on the lattice. 
Thus, a fundamentally new approach to the long-standing problem of structure functions 
at very small xbj seems possible. We give arguments that the limit xbj corresponds to 
a critical point of the effective field theory where the correlation length becomes infinite 
in one direction. 
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1 Introduction 



The behaviour of the structure functions of deep inelastic lepton-nucleon scattering (DIS) 
at small values of the Bjorken variable xbj has been a topic of intense experimental and 
theoretical research in recent years. The strong rise of the structure functions for xbj — * 
observed at HERA opened the way for a more detailed experimental study of the 
small-XBj limit of QCD. Today, many theoretical models try to explain the behaviour of 
the structure functions in this limit. 

In a very pragmatic perturbative approach, the well-known DGLAP dynamics of the 
evolution of structure functions is taken seriously in the small-xej region (cf., for 
example, the parametrizations of [Q). Given a sufficiently small starting scale of the 
perturbative evolution, the small-xej rise of structure functions can even be described 
on the basis of valence-like parton distributions (see ^ and refs. therein). In this case, 
all of the small-XBj rise is produced by perturbative QCD. 

If, in a more modest approach, a constant behaviour or a soft growth for Xbj — > 
is ascribed to the input distributions, the DGLAP evolution allows for a good fit to the 
data even in the case of a larger starting scale. This observation underlies successful 
analyses based on the idea of double asymptotic scaling 0]. It is also at the heart of 
a recent combined analysis of diffractive and inclusive structure functions P], where a 
universal logarithmic growth of the input distributions was assumed. 

An obvious problem of the above, purely DGLAP-based approaches is the presence 
of large ln(l/a;Bj) corrections, which can spoil the usefulness of the perturbation series as 
an asymptotic expansion. Using the BFKL resummation method (see [Q for reviews 
and further references), all leading ln(l/xBj) contributions can be included in the struc- 
ture function analysis (see, e.g., |^). It is, at present, not completely clear whether the 
predicted very strong power-like growth is borne out by HERA data. Also, no final con- 
clusions can yet be made concerning the effect of the recently obtained next-to-leading 
order results in the BFKL framework (see [|10| and refs. therein) on structure function 
analyses. Different methods of dealing with the apparent large size of the corrections 
have been suggested [pH) . 

Unfortunately, because of the problem of infrared diffusion [jl2|, it is not clear that 
the perturbative resummation of ln(l/a;Bj) terms at any given order in can reveal the 
asymptotic behaviour. One possibility to go beyond the above resummation schemes is 



the construction of a high-energy effective action for reggeons (see [13| and refs. therein). 
Another possibility is the derivation of a small-XBj evolution equation in the framework 
of the semiclassical approach. This method goes back to the treatment of hadron-hadron 
scattering in the functional integral approach using the eikonal approximation introduced 
T^. For recent results concerning the energy dependence of cross sections in this 



m 



framework see, e.g., [[T5|] . 



A very different point of view is advertised in |TB], where, in addition to the familiar 
soft pomeron, a phenomenological hard pomeron is introduced to describe small-XBj 



structure functions. A further analysis motivated by the successful phenomenology 



of |jT6|, concludes that the appearance of new small-XBj singularities due to evolution 
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is inconsistent with the analyticity principle governing Regge theory and underlying the 



approach of [If:]. Note also the successful analysis of various experimental cross sections 



reported in [jT8[, which combines ideas of fT^ and [|16 



In spite of the large amount of work invested and the good description of data 
achieved in most of the above approaches, our fundamental understanding of the small- 
xbj asymptotics in QCD is still not satisfactory. Therefore, we propose yet another 
method to achieve this goal. In this article, we will derive a connection between the 
small-XBj behaviour of the structure functions of DIS and the behaviour of correlation 
functions in a certain effective field theory, which we consider both in Minkowski and Eu- 
clidean space-time. The analytic continuation from Minkowski to Euclidean space-time 



to be discussed below is similar in spirit to the one described in ||T9[ for the high-energy 
hadron-hadron scattering amplitudes. However, the details are quite different and we 
will comment on the connection between the two approaches. 

For simplicity, we study in this article a scalar model field theory where scalar cur- 
rents replace the electromagnetic currents of real life, i.e., QCD. However, it will become 
clear that our methods allow a straightforward generalization to QCD. 

Our article is organized as follows. In Sect. |^ we define the model and discuss the 
properties of our scalar analogue of the virtual Compton amplitude. In Sect. ^ we in- 
troduce a class of effective Hamiltonians depending on a parameter r. We also make a 
rotation from Minkowski to Euclidean space-time and show that the high-Q^, xbj 
limit corresponds to r — 0. In Sect. ^ the amplitude is written as a functional integral 
with r-dependent action. In Sect. § the formalism is applied to a simple model of free 
fields. We show that - at least in the free-field limit - we get a diverging correlation 
length in our effective theory for r — > 0, i.e., a critical phenomenon. Finally, two alter- 
native continuations from Minkowski to Euclidean space-time are outlined in Sect. |^. 
Section |^ contains our conclusions. 



2 The model and the scalar analogue of the virtual 
Compton amplitude 

We consider a model with one real scalar field and the classical Lagrangian 

£{x) = - d,<P{x) dy{x) - — - ^ , (2.1) 

where m and A are the unrenormalized mass and coupling parameter, respectively. The 
momentum canonically conjugate to (pix) is 

djC 

vr(x) = = 0(a;) . (2.2) 
The Hamiltonian and the third component of the momentum operator read 

H = J^^^^^^^^ d^xi^^Uixy + l{VHx)f + "^^xY + A$(x)^| , (2.3) 
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In Eqs. ( p.3|) and (|2.4|) no vacuum expectation values are subtracted. Note that we 
use lower- and upper-case letters for the classical field variables and the corresponding 
operators respectively. 

We assume that the theory has physical states with reasonable spectrumQ etc. We 
will call 'proton' the particle of lowest mass M. We introduce now an 'electromagnetic' 
coupling of our field to a scalar 'photon' with field A{x): 

C\x) = -e J{x) A{x) with J(x) = . (2.5) 

This scalar current is our analogue of the electromagnetic current of QCD. In the fol- 
lowing we will study DIS, i.e., the cross section for the absorption of the scalar 'photon' 
on the scalar 'proton' (Fig. |l|). 




^ ^^P'^ 

P 

Figure [T|: Diagram for the absorption of a scalar photon on the scalar proton. 
We define the structure function for this reaction as 

Wiu^Q') = E k2nYS^'\p'~P-q) (P(p)|J(0)|X(p')) {Xip')\.m\Pip)) , (2.6) 
X ^ 

where 

= -g^ > and z/ = j9g/M>0. (2.7) 

All the notation is standard. The conventions for the metric, normalization of states etc. 
follow |20i. 



The object to study from a theoretical point of view is the virtual Compton ampli- 
tude. To be precise, we will study the amplitude corresponding to the retarded commu- 
tator (see, e.g. 



TAv.Q^) = —J d'xe^'^^eixo) (P(p)|[J(x), J(0)]|P(p)) . (2.8) 

Here and in the remainder of the paper, we always assume that the connected part of 
the relevant matrix element is taken, (■■■) = (■■ ■)c- It is easy to see that 

ImT,(z/ + le, Q") = sgn(z/) W{\u\,Q^) (2.9) 

for u real and e —>■ 0+. 



^ Of course we know that the theory defined by Eq. (2.1) is trivial if the cutoff goes to infinity in 



four dimensions. The purist may always think of a theory with a finite cutoff. 



4 



We consider now Eq. (|2.8|) in the proton rest frame, where we can set 



p={M,0) , q=(u,es^u^ + Q^) (2.10) 

with Cj {j = 1,2,3) the spatial cartesian unit vectors. Rotational invariance gives 

(P(p)|J(f,xO)J(0)|P(p)) = (P(p)|J(±|f|e3,x°)J(0)|P(p)). (2.11) 

Inserting this in Eq. ( |2.8| ) and using Eq. ( [^.1(J| ) and the fact that the commutator in 
Eq. (p.8|) vanishes outside the light cone, we find 



Tr{u, Q') = - r^^—^ J dx' J dx'x' ^i^^u-^.^V^' {P{p)[J{x%, X°), J(0)]|P(p)) . 

(2.12) 

In Eq. ( |2.12| ) we have a convenient starting point for the analytic continuation of 
Tr{h',Q'^) into the upper half of the complex u plane, keeping always fixed. In fact. 



it is easy to see that the integral in Eq. (2.12) defines an analytic function for Imz/ > 



(cf. Appendix^. Using standard methods, one finds then that Tr(z/, Q^) can be extended 
to a function analytic in the whole u plane except for poles and cuts on the real axis for 



> gV2M (cf. Fig. i). 



/ 


\ Imv 


- Q^/2M 


Q^/2M 



Rev 



Figure The analyticity structure of the virtual Compton amplitude Tr{i',Q'^) in the 
I' plane. Singularities can only occur for real z/, |z/| > Q'^/2M. The integral in Eq. ( |2.12| ) 
is a representation of Tr(z/, Q^) for Imz/ > 0. 

We are interested in the behaviour of Tr(z/, Q^) for ^ oo on the real axis. Instead of 
studying this directly, we will study Tr(z/, Q^) on the imaginary axis, i.e., for u = ir], r] ^ 
+00. With standard assumptions for the high-momentum behaviour of the amplitudes 
of QFT (cf. Appendix 0), we can apply the Phragmen-Lindelof theorem (see, e.g., [0), 



which ensures that the asymptotic behaviour on the real and imaginary axes are related 
by analytic continuation. Our object of study is thus: 

oo x° 

Tr{tr],Q') = /_ / J cia;Ve-^"''+^'v/^(P(p)|[J(a:3e3,x°), J(0)]|P(p)). 



-xO 



(2.13) 

It turns out to be convenient to generalize this amplitude slightly. Let fi with Im /i > 
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be a complex parameter with dimension of mass and define: 



(2.14) 

We have 

r,(^r/,Q2) =ReT+(«r/,Q2^0), (2.15) 

which is obvious for rj > Q and requires a httle thought for < r/ < Q (cf. Appendix ^). 

In the next section we will represent T+(ir7, Q^, /i) for Re/i < — M as a correlation 
function in an effective Euclidean field theory. The amplitude Tr{iri, Q^) is then obtained 
from Eq. (|2.15| ) after analytic continuation of T+ to yU = 0. We will also write T+{iri, Q^, 0) 
as a correlation function in the corresponding effective theory in Minkowski space-time. 



3 From Minkowski to Euclidean space 



Let us introduce new coordinates in Eq. ( |2.14|) : 

x° = ^ cosh X , = ^ sinh x 

with 

< ^ < CX3 , —00 < X < 00 . 



(3.1) 



(3.2) 



In this section we always consider t] > Q. Then we insert a sum over a complete set of 
states and perform the ^ integration: 



T+{i7],Q'^,fi) = J2 _ Q2 J^^ dx sinh x 



(3.3) 



xe 



( r?+i (p' ° - Af -^t ) ) cosh X - ( ^Jrp^^^Viv' ^ ) sinh x 



(P(p)|J(0)|X(p'))(X(p')|J(0)|P(p)). 



The convergence of the ^ integration in Eq. ( p.3| ) is guaranteed since 



ri cosh X^ yff' ~ sinh x > . 



and we get 



T+{iri,Q^,fj.) = / 2_Q l J_^^^ ^^^^^ 



X 2 I (^?7 + i(p'^ — M — yu)^ cosh x — (y\frf^-Q^ + ) sinh x 

^{p{p)\jmnp'))W)\jmp{p)) 



(3.4) 
(3.5) 

-3 



Since p'^ > \p'^ \ , we have also 



{p^ — M — Re /i) cosh x ~ p'^ sinh x > 



(3.6) 
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for Re/i < —M. This allows us to write Eq. ( p.5| ) again as an integral over a parameter ^: 

-2 



T+{ir],Q'^,fi) = J2 



dxsinhx / d^^ 
Jo 



(3.7) 



X e 



x(P(p)|J(0)|X(p'))(X(p')|J(0)|P(p)) 
-^^==,JjXsn.hxl die 



X e 



{fi+ir}) coshx— Tj'^—Q'^ sinhx 



= (P(p)l eC(Hcoshx-P3 sinhx) J(0) e-?(/^coshx-P3 sinhx) J(0)|P(j9)) . (3.8) 



where 



Effectively we have shifted the ^-integration in Eq. p.3| ) from the real to the imag- 
inary axis, i — > —2^. In Eq. ( p.8|) H and P^ denote the Hamilton operator and the 
third component of the momentum operator of the original theory, Eqs. (|2.3|) and ( p.4|) , 
respectively. 

The matrix element e can now be interpreted as one of a Euclidean field theory. 
For this we set 



and 



We have 



2/° = ^ cosh X , = ^ sinh x 

y+ + y- y" 



y+>0, y->0 , < r < 2 . 
Furthermore we define an effective, r-dependent Hamiltonian: 

H,ff{r)=H-{l-r)P\ 



(3.9) 

(3.10) 
(3.11) 

(3.12) 



With this we can write 

Me ^ ME{yo, r) = (P(p)l e^"^^//W J(0) e-^°^^//« J(0)|P(p)) (3.13) 

and 

T+(^r?,g^/i) = ;^^=^= / dy+y^ d2/_(y+-y_)e^^+(^+^^-)+^^-(^+^^+)A^^(2/°,r) 

(3.14) 

V(/)' f^^(l - r) e?'''[A'+*(i~i)'?-+^i^+l A<^(?/0, r) . 



_2 ^oo 
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We see from Eq. ( 3.1CI|) that for rj — > cxo, fixed, we have rj^ ~ 21], ~ Q^/ {2r]). Then 



the oscillating term exp{iy^rri^/2) restricts the r-integration to 

r<^. (3.15) 

vy 

Thus, except for the small- range, 0<y^<l/r], the behaviour of the matrix element 
A^^(?/o, r) for r will be essential. Going from r = 1 to r = in Eq. (|3.12D corre- 



sponds, of course, to going from the ordinary Hamiltonian H to the light-cone Hamilto- 
nian H — P^. In conventional light-cone calculations the theory is quantized directly on a 
light-like subspace (for reviews and further refs. see p3[)- However, we always stay away 



from the light cone by having r ^ 0. For studies concerning the possibility to approach 
the light cone continuously in the Hamiltonian method, we refer to [Q. In Sect. ^ be- 
low we will use the Lagrangian method and for us the light-like limit concerns only the 
endpoint r = of the r integration in Eq. ( p. 141 ). 



Whereas for r G (0, 2), the effective Hamiltonian Hf.ff{r) has an energy gap, this gap 
vanishes for r = and 2. We have 



mm{{X\H,ff{r)\X) - Eo} = ^r{2 - r) M , (3.16) 

where Eq is the vacuum energy: 

Eo = m,ffir)\0) = {0\H\0), (3.17) 

and the minimum is taken over all normalized states |X) orthogonal to |0). Equa- 
tion ( p. 161 ) is already indicative of a large correlation length and a critical phenomenon 



for r = and we will give more arguments for this in Sect. 

Instead of representing T^{iri, Q^, fi) as an integral over a Euclidean matrix element 



A^^(?/°,r) in Eq. ( |3.14| ), we can also represent it as an integral over the corresponding 
Minkowskian matrix element M.E{y^if^) = ■M.M{—W^,f')'- 



oo 2 

2z 



T+itV,Q',f^) = ^J^T^ I dx\x'f I rfr(l-r)e-^°H'^+(^-i)^-+^^''+]^M(a^°,^)- 







(3.18) 

with 

X^,,(xO,r) = (P(p)|e+'"''^«//(^) J(0)e-'"°^=^^(^) J(0)|P(j9)). (3.19) 
In Eq. ( p.l8| ) there is no problem when setting /i = 0. 



4 Functional integral representation 

In this section we will rewrite the matrix element of Eq. (|3.13|) in terms of a path integral, 
using standard procedures as in the conventional case -ffe// = i.e., r = 1. We assume 
for simplicity that the states \P{p)) have the quantum numbers of the $ field, i.e., that 
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the $ field is an interpolating field for these states. Then, according to the basic principles 
of the LSZ reduction formalism, we have to construct the operators 

A(p, x^)=i I ^xe'P'^ tAt "^(a;) • (4-1) 

JxO=const. 

We have then, in the weak sense, 

lim ^At(p,x°)|0) = \P{p)Z') = \P{p)) , (4.2) 

where Z is the wave function renormalization constant. 

The Euchdean version of Eq. ( [4.2|) for the proton with zero three-momentum reads: 

1= hm (e-'^'AUom) e(*^+^«)* = \P(p)) , (4.3) 



where we set A{0) = A{p,0) with p = (M, 0). In terms of the field operator and its 
conjugate canonical momentum, we have 

A{0) = [ d^x{m{x) + M^{x)) , (4.4) 
A^O) = [ d^x{~m{x) + M^{x)) . (4.5) 



Before rewriting the matrix element of Eq. ( p.l3| ) in terms of a path integral, let us 
recall the corresponding procedure in the conventional case -f^e// = H, i.e., r = 1. Using 
Eq. (Of) we have 



MEir,l) = - lim e^-f-')i^+^o) /0U(0) e-^'^"")^ J(0) e-("-°)^ J(0) e-(°-"')^ At(0)|0) . 
Z 

(4.6) 

Now the standard procedure is to introduce alternating intermediate states which are 
eigenstates of the field operators 11 and $ with definite classical field eigenvalues vr and 
(/) on a sufficiently dense grid in Euclidean time. The vr integrations can then be carried 
out explicitly, leaving one with a product of integrations. In the limit where the time 
slices become arbitrarily thin, this defines the Euclidean path integral. 

In realizing this procedure, all operators in Eq. ( f4.6| ) are replaced by functionals of 
vr and (j). For H and J, the relevant expressions are given in Eqs. ( p.3|) and ( p.5|) , for the 
annihilation and creation operator A, A'^ in Eqs. ( [4. 41 ) and ( (4.51 ). 

The standard procedure outlined above gives the result 

Me{tA) = ^ lim e^^f-^^^'' f D<Pa{Tf)j{T)j{0)a\n)e-I''''^- , (4.7) 
Z J 

Z = J D^e-/"^*^^^ , (4.8) 
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with the Euchdean Lagrangian 

£e = -i7i{x)4){x) +n (4.9) 
= -nx{x)^{x) + -7i{xf + - (V0(a;)) + —<p{xf + -<P{xf . 



Here the canonical momentum vr entering Eq. ( |4.7|) via a, and Ce has to be understood 
as a function of 0. In the Euchdean theory, the relation between cf) and vr is defined by 

rll-l 

i<P= — = 7r. (4.10) 

(777 



In this way, Eq. (|4.71 ) and the expressions for a and follow directly from Eq. ( [4.6|) 
and the explicit form of the Hamiltonian density Ti. 

The result at r 7^ 1 can be obtained by following the standard procedure for convert- 
ing a matrix element into a functional integral expression outlined above. It is easy to 
guess the correct answer by observing that the effective Hamiltonian density is defined 
by the substitution 

n — > Heff = n + (1 - r)7ids<l) . (4.11) 
Thus, the expression for tt in terms of (f) following from Eq. ( [4.101 ) now reads 

^^=^^ = 7i+{l-r)^s<p. (4.12) 

This is to be used when expressing a, and the new effective Lagrangian C-E,eff in 
terms of (f) in Eq. (|4.7|) . 

We summarize the results of this calculation, some details of which are given in 
Appendix ^ the Euclidean matrix element of Eq. ( p.l3|) can be calculated as 

ME{T,r) = ^ ^lirn^ e^rf-n)M^-i J D(j) a{Tf)j{T)j{0)a\Ti) e' I '^'''^''■^ff , (4.13) 

Z = J D(j)e-I'^''"^^'^ff , (4.14) 



Tf — >- + 00 



with 



1 1 Z \ 

^E,eff = ^(5o0(x) + z(l - r) ds^y + -(V0(x))' + ^<P{xf + -0(x)^ (4.15) 



and 

a{Tf) = f d^x(-do(pix)-i{l-r)d3(p{x) + M(l){x)), (4.16) 
a^{Ti) = [ Sx [do(j){x) + t{l - r) 930(x) + M0(x)) , (4.17) 



j(r) = 02(O,r). (4.18) 

Note that in the Euclidean path integral, Eq. ( |4.13| ), the functionals a(r) and a^(r) are 
not complex conjugate to each other and C-E^eff is not a real function for r 7^ 1. 
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The matrix element Eq. ( 4.13 ) is related to the high-energy limit of the structure 
function W via Eqs. ( p.9|) , ( p.l5|) and (|3.14|) . The above Euclidean functional integral 
expression for A^^; is one of the central results of this paper. We propose to evaluate it 
using genuinely non-perturbative methods, e.g., lattice Monte Carlo simulations. 

Similarly, the Minkowskian matrix element of Eq. (|3.19| ) can be expressed in terms 
of an appropriate Minkowskian functional integral (cf. Appendix 0). 



5 Free-field example 

To illustrate the meaning of the formalism developed, let us consider the simple example 
of free fields, i.e., the case A = 0. 

Using Eqs. (|4.6|) , ( [4.16| ) and ([4.17|) , the matrix element is expressed as 



ME(T,r) 



where 



- hm e^-f-^)^' [£^jd'xJ-^ + 
Z n^-^ J \ drr 



Xi = {Xi,Ti) , Xf = {Xf,Tf) , X={0,t) 



(5.1) 
(5.2) 



and Gfi^jjff, is the connected Green function of two fields and two currents j illustrated 
in Fig. 0. The terms with x^ derivative disappear because of the x^ integration. 



j(x) 





j(0) 



Figure |^: The Euclidean Green function required for the calculation of Ale . 

In the free-field case. Fig. ^ contains only the two diagrams shown in Fig. ^, where 
the propagator following from the free part of the Euclidean Lagrangian, Eq. ([4.15| ), has 
to be used. It reads 



G{x) 



Akx 



m 



'X 



Ki ( mJx'^A^ ^x 



(5.3) 



where 



f 1 

10 

l-(l-r)2 i{l-r) 

V i(l -r) 1 



with X 



f x^ \ 

x^ 
x^ 
V r ) 



( k^\ 
k2 
k3 
\h J 



(5.4) 



Ki is the modified Bessel function, and m = M m our free-field model. It is easy to see 
that x'^A~^x ~ 2rr^ for x = (0, r) and r 0. According to Eq. {f).3\j, this means that 
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Figure |^: The two diagrams contributing to the Euchdean Green function Gfpjjrp in the 
free-field case. 



the correlation length in the Euclidean time direction diverges and we are dealing with 
a critical phenomenon. 

To study in general the behaviour of G{x) for a; ^ oo we set 



X = Rx , R = \x\ = y xf + X2 + xl + t"^ (5.5) 
and consider the limit i? — 00 for fixed unit vector x. We get from Eq. ( |5.3D : 

G{x) {m^R'x^A;'x) i? ^ 00 . (5.6) 

This 



1 fW 

Thus the decay of |G'(x)| as i? — > 00 is governed by exp —MR Re^x'^A~^x 
suggests to define a correlation length depending on x as 



e(x) 



■1 -1 



MReJxTA-^x 



We have, in the Euclidean time direction, 



e(o, 1) 



Mx r(2 - r] 



(5.7) 



(5. 



and, in the space directions, 

e(e„0) = M-^ (5.9) 

Thus for r — > the correlation in this free-field model becomes of long range only in the 
time direction. There ^(0, 1) diverges as I/a/t, i.e., in a manner characteristic of mean- 
field theories for critical phenomena in statistical mechanics, identifying r with T — Tc, 
the deviation of the temperature from the critical value. 

We also note the behaviour of G{x), Eq. (|5.3|), for — > 0: 



G{x) 



— (r^xA^^ 



-1 



for 



R^O. 



In the Euclidean time direction this can be written as 

^(0»^^f^V for .«^(0,1) 



(5.10) 



(5.11) 



This means that for times much smaller than the correlation length the Green function 
in the free theory has a simple power behaviour, again similar to what is found in the 
study of critical phenomena. 
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Now we describe the calculation of in Eq. ( |3.14|) in the free-field case. Changing 
the integration variables in Eq. ( p.l4| ) from y± to and y^, using the expression of 
Eq. (|0| ) for and calculating the Green function G^jj^ according to Fig. ||, the 

amplitude T+ reads 



T+{iri,Q'^,fx) 



3 „,3 ^|[(i/0+y3)(M+ir,_)+(yO-3/3)(M+ir,+ )] 



X 



M 



Mv/(y0)2-(y3)2 . (5.12) 



Applying the identity 



K, [M^iyoy - (1/3)2 



d 



Ko MJ{yO)2 - (^3)2 



integrating by parts, and returning to the variables y±^ one finds 

j roo roo 

T+(ir/,Q^/i) = ^ dy+ ^-,[y+{^^+iv-)+y-{^^+^v+)] 

477^ JO Jo 



(5.13) 



gM(s,+ +j/_)/2 _^ g-M(s,+ +y_)/2 



(5.14) 



For sufficiently negative Reu, the integrals can be evaluated using the relations 6.631.3, 



9.224 and 6.224.1 of [251. The resuh is 



T+{zV, Q\ /i) = {/(M) + /(-M)} 



with 
/(M) 



1 



/i2 + 2M/i - g2 + 2iri{M+ij) 



In 



_{(M+/i)2-g2 + 2^r/(M+/x)} 



(5.15) 



. (5.16) 



When this expression is analytically continued to /i = 0, a careful treatment of the 
imaginary part of the logarithms is essential. The singularities of T+(i77, Q^, /i) in the 
complex fi plane are shown in Fig. ^ The continuation in fi from Re < —M to /i = 
has to be done along the real axis staying above the left hand singularities. 

From Eq. ( p.l5|) one finds the behaviour of the retarded amplitude on the imaginary 
axis in the u plane: 



Tr{ir],Q^) = ReT+{ir],Q^,0) 



+ 



27V \Q^- 2ir]M + 2ir]M 



(5.17) 



The corresponding behaviour on the real axis follows from the analytic continuation 
17] ^ V for V real: 



1 



1 



+ 



27r \Q2 _2zyM-ie Q^ + 2uM + ie 



(5.1^ 
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M 



M-ir\ 
M-ir[_ X 

M-ir\^ X 



M 

h- 

M-ir\_ Xr- 
M-ir[_ X 

M-/rj+ X 



-> Re|i 



Figure ^ The singularities of T+(ir7, Q^, /i) in the complex fi plane. There are poles 
at ±M — ifj^, ±M — ifj_ with fj± = t] ± y/rj'^ — Q'^ — M'^ and logarithmic branch cuts 
starting at ±M — irj^ and ±M — iri_. Here we assume yf > Q"^ + M^. 



This is in agreement with what one would have found by simply adding the two dia- 
grams for the forward scattering of a scalar proton and a scalar photon in the original 
Minkowskian theory. 

We will close this section with some speculative remarks. Looking at Eq. ( |3.14| ) we 
see that for ^ oo the typical integration range over M-Eiy^^f) is 

< r < ^ , (5.19) 



Q<y'<^^yl. (5.20) 



2r] 

Let us insert the mean value of ip of Eq. ( p.l4| ) into Eq. ( |5.19| ) to get as typical value 
for r: 

f=^. (5.21) 

Then we have two relevant scales for the integration in Eq. (|3.14| ): and the length 
scale of the decay of M.E[y^ Also in the general case of the theory with interactions 
this should be governed by a correlation length ^{r) similar to ^*^°^(f) = ^(0, 1), Eq. ( p.8| ), 
in the free-field case. Clearly the behaviour of the amplitude T4. will depend crucially on 
whether is smaller or bigger than ^(f). 

In the free- field case we have 

< e^°)(f) for Q2>8^2_ ^5 22) 

Here for large and r/ — > cxd the integration in Eq. ( |3.14|) probes only the region of 
M.Eiy^^'r) where is smaller than the correlation length and thus where Eq. ( |5.11| ) 
applies. We can speculate that in the theory with interaction we should again be able to 
distinguish two regimes, 

yl. i m , (5.23) 
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where, from our experience with critical phenomena, we would expect ^(r) to satisfy a 
scaling law for f —>■ 0, 



c = const. 



p > 0, 



(5.24) 



but not with the mean field exponent p = 1 /2 as in Eq. ( |5.8|) . 

For < p < 1 the dividing line between the two regimes of Eq. (|5.23|) would then be 
roughly: 

^ = — I2M j ^'-2') 
and it is tempting to speculate that this could correspond to the dividing line between 
hard [y'^ < ^(r)] and soft [y^ > ^(r)] pomeron regimes. But much more work is needed 
before such a conjecture can be substantiated. 

Let us, nevertheless, assume for the moment that such a conjecture is true and apply 
Eq. (|5.25 ) also in the physical region, replacing rj by u. Then we should find a dividing 
line between hard and soft pomeron effects in the structure functions in the xbj versus 

plane: 



(5.26) 



This is illustrated in Fig. ^. Thus, taking the limit xbj — * at fixed > we would 
finally always go from the hard to the soft regime for < p < 1/2 and vice versa for 
1/2 < p < 1. 



^ Q' 



< p < 1/2 



hard 



soft 



-> X 



Bj 



1/2 < P < 1 



hard 



soft 



Bj 



Figure |^: Sketch of the regions of soft and hard pomeron effects in the xbj-Q^ plane 
using Eq. (|5.26| ). 



6 Alternative analytic continuations to Euclidean 
space 

It is clear that the form of the analytic continuation from the Minkowski to the Euclidean 
theory discussed in the previous sections is by no means unique. Other possible forms of 
analytic continuations from Minkowski to Euclidean space can be derived: in this section 
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we discuss two of these alternative Euclidean formulations and compare them with the 
one of Sect. |^. 

Working in the rest frame of the proton, consider the quantity 

= ^ [ d'xe{x')d'^%P{p)\J{x)JmP{p)) (6.1) 

Ztx J 

as an analytic function of with the spatial part q fixed and real. Using standard 
methods, one finds 



nq',q) = - E(2vr)^^(^Hg - p') n ^ ,\ _ J (P(P)\J(^)\^(P'))\" ' (6-2) 
X q -\- ivi p 



which reveals the cut along the positive real axis. The structure function of Eq. ( p.6|) is 
related to T by 

lm[Tiq' + ie,q)]=Wiu,Q^) , (6.3) 

where e —>■ 0+. Now, consider T{q^,q) on the positive imaginary axis, i.e., for g° = irj 
with 1] > 0. Using the identity 



ir] + M - p'O 



-- „'0 



dx4e('''+*^-^' (6.4) 



(where the problematic lowest value of p'^, p'^ = M, can be treated by giving rj a small 
positive imaginary part), the following expression for T is obtained from Eq. 



1 f f + OO ^ 

T{^V,q) = — / / dx,e<^^'-''-^\P{p)\JE{x,x,)JEmP{p)) ■ (6.5) 

ZTT J Jo 

Here Je{x,X4) is the current of the conventional Euclidean theory: 

Je{x, X4) = e^^V(0, f)e-^^* . (6.6) 

Using the methods of Sect. ^ [cf. Eqs. ( p.llj ) and ( p.l2| )], invariance under spatial rotations 
can be employed to further simplify the above expression for T. 

To summarize, the amplitude T on the real axis, which determines the physical 
structure function, can be obtained from the one on the imaginary axis by analytic 
continuation in (g° = irj) — »• (g° = ly). According to Eq. ( |6.5| ), the latter one can be 
calculated from a conventional Euclidean matrix element. 

This approach is simpler than the one discussed in the previous sections in that the 
Euclidean theory is based on the usual Hamiltonian H and not on the more complicated 
effective Hamiltonian of Eq. (|3.12|) . However, in order to derive the large-z/ behaviour 
of the structure function W{i>, Q"^) at fixed Q^, one needs the analytic continuation of 
T{q^,q) in for every value of q. Only at this point can one study the asymptotic limit 
u — > 00, which involves both q^ —>■ 00 and |g | — > 00 in the physical region. This analytic 
continuation in g° is expected to be more troublesome than the /x continuation of Sect. ||, 
which is decoupled from the high-energy limit. 
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We consider now a third possible form of analytic continuation from Minkowski to 
Euclidean space, which is close in spirit to ||19|. The basic object is still the virtual 
Compton amplitude Eq. ( p.8|) , which we write now as: 

Triu,Q') = l{f^iq)+fl{-qn) (6.7) 



where 



f+{q) = - l(fxe{xy{{x''f-{x^f)e"^^P{p)\J{x)J{<S)\P{p)). (6J 



TT 

Here we work in the rest system of the proton and choose now q in the form 



q={u,e^^p^ + Q\Q), (6.9) 

with e± a unit vector in the {q^iq^) plane. With the parametrization x = 
(^ coshx, x±, ^ sinh x), we find 

TT J J-oo Jo 

(6.10) 

We will show how to compute this quantity and its u derivatives at z/ = using an 
analytic continuation from Minkowski to Euclidean space. First of all, we write 



T+(g)|,=o = / d'x^e-"^--'- / rfxA«(x,f±) , (6.11) 
where q± = Qe± and the function A^^^ is defined by 

A^'\x,x^) = - / d^^{P{p)\J{x)JmPip)) 

TT JO 

= -T.{Pip)\Jio,x^Mx{p')){x{p')\jmp{p)) 

^ X 

X 2 > (6-12) 

[{p'^ — M) cosh X — p'^ sinh x ~ i^] 

where finally the limit e — > 0+ is taken. The singularity structure in the complex x plane 
is such that the function can be continued to the imaginary axis, x ~^ where, however, 
the 6 range is restricted by 6* G (— 7r/2, 7r/2). In this region, A^^"* is related to a Euclidean 
matrix element, 



A^'\te,x^) = — / d^^{p{p)\MxE)JEmp{p)) 

TT JO 

= -E(^(p)I'^(o,^±,o)|x(pO)(x(p')|J(o)|p(p)) 



^ X 



X 2 , (6.13) 

[(p'O — M) cos 9 — ip'^ sin 9 — ie] 

where xe = (x±,^sm9,^cos9) is a Euclidean four- vector and the Euchdean current 
Je{xe) is defined by Eq. (|6.6|). 
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More generally, the expression for the n-th derivative of T+ with respect to p, eval- 
uated at z/ = 0, reads 



" f r+oo 

a,"T+(g)|,=o = E / d'x^e-"'--^-Cr{q±-x^,Q') / cix(coshx)'A('+i)(x, f^) , (6.14) 
with explicitly calculable coefficients and functions A^'^ defined by 



T.iPiP) I -^(0, 0) \X{p')){X{p') I J(0) |P(p)) 



X ^ m • (6.15) 

[{p'^ — M) cosh X — p'^ sinh x — i^] 



As in the case / = 1, the values of this function for real x can be obtained from 

A^'\^9,x^) = r die{P{p)\JE{xE)JEmP{p)) 

TT Jo 

E (P(P) I ^(0, 0) \X{p')){X{p') I miPip)) 



TT Jo 



^ X 



X . . „ . . . ,^+1 • (6-16) 



1 

[(p'O — M) cos 6* — ip'^ sin 6* — ie] 
by the analytic continuation i6 x- 

To summarize, one first computes A'^^\i9, x±) for / = 1, . . . , n+1 and 9 E (— 7r/2, n/2) 
in the Euclidean field theory (for example, on the lattice). Then, the value of d^T+{q)\^=o 
can be obtained from Eq. (|6.14| ) using the analytic continuations of the A^''^ to the 
real axis. Thus, using Eq. ( |6.7| ) and a Taylor expansion of T+ around 1^ = 0, the be- 
haviour of Tr(z/, Q"^) at small u (corresponding to large xbj) can be studied. The values of 
d'^Triv, (5^)|i/=o can then be related to integrals of the structure function Wiv, Q"^) using 
dispersion relations in v. Assuming for simplicity no subtractions this reads: 



1 r'^ 

Tr{u,Q^) = - du'W{iy',Q^' 

71 Jq2/2M 



1 1 



u — u — u' + u + le 



(6.17) 



ril /-oo W(u' O"^] 

d^Triu,Q^)\^.o = -l^^^^^du'^^ (n = 0,l,2,...). (6.18) 

Again, the advantage with respect to the method of Sect. |]is the use of the conventional 
Hamiltonian in the Euclidean theory. The disadvantage is the limited sensitivity to the 
region of small Xbj. 



7 Conclusions 

In this paper, a fundamentally new approach to the long-standing problem of small-xej 
structure functions in DIS has been developed. It is based on the well-known relation 
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of the small-XBj limit of DIS and the high-energy hmit of forward virtual Compton 
scattering. Instead of taking the energy of the forward virtual Compton amplitude to 
infinity along the real axis, we propose to consider the limit of large imaginary energy. 
According to the theorems of Phragmen and Lindelof type, the asymptotic behaviour in 
both directions is related by analytic continuation. 

A slight generalization of the above virtual Compton amplitude with imaginary en- 
ergy, which contains an additional mass variable /x for convergence, can be written as an 
integral over a matrix element in an effective Euclidean field theory. Using standard meth- 
ods, this matrix element is expressed in terms of a Euclidean functional integral with a 
simple effective Lagrangian, which is given explicitly. The essential proposal of our paper 
is to attempt an evaluation of this functional integral using genuinely non-perturbative 
methods, e.g., lattice Monte-Carlo simulations. Finally, an analytic continuation to yU = 
and to the real axis in the complex energy plane should allow the extraction of the desired 
asymptotic small-xej limit of DIS structure functions. 

To illustrate how our approach works, we have explicitly performed all necessary 
steps in a simple model with free scalar fields. In this case, the Euclidean path integral is 
given by two tree-level Feynman diagrams and the subsequent analytic continuations in /i 
and the energy variable can be carried out explicitly. As expected, the final result for the 
structure function is in agreement with a direct diagrammatic calculation in Minkowski 
space. 

We have shown that in the free-field case the high-energy limit of the virtual Compton 
amplitude is governed by the behaviour of the effective Euclidean theory near the critical 
point r = 0, where r is our parameter defined in Eq. ( p. 10 ). We have then assumed that 



also in the effective theory with interaction a critical point occurs at r = and discussed 
the possible consequences of this for the hard and soft pomeron regimes of the structure 
function. 

We do not claim that we have found the optimal method for an analytic continuation 
of high-energy amplitudes. In fact, in Sect. || of this paper we introduce two alternative 
possibilities. They have the advantage that the Euclidean field theory obtained is based 
on the conventional Hamiltonian as opposed to the effective Hamiltonian required in the 
original approach. However, in both cases it appears to be more difficult to recover the 
small-XBj limit of structure functions from the final result of the calculation. 

A lot of work remains to be done before phenomenologically relevant information 
can be extracted from the approach suggested. Note first that all of the discussion in the 
present paper is based on a model with a scalar 'photon' coupled to scalar partons. This 
has to be extended to the realistic case of a vector photon and QCD. However, we do 
not expect any fundamental problems with this generalization. 

Furthermore, the feasibility of a lattice calculation or any alternative non-pertur- 
bative treatment of our Euchdean path integral expression has to be evaluated. Note 
that our effective Lagrangian contains an imaginary part, which, even though it does not 
raise any fundamental problems, may lead to technical difficulties on the lattice. Note 
also that the subsequent analytic continuation back to the physical region could prove 
highly non-trivial. 
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On a more fundamental level, one has to ask whether better methods for the treat- 
ment of the small-^Bj limit of DIS in a Euclidean theory exist. We have attempted to 
address this question in our cursory discussion of alternative methods of analytic con- 
tinuation. However, we are not able to give a general answer at present. 

Given the above reservations, the importance of our results lies in their potential to 
relate small-xej cross sections in DIS to quantities derivable from an effective Euclidean 
field theory. One can then hope to understand the limit xbj — > in terms of the critical 
behaviour in this Euclidean theory using all the tools which have been developed in sta- 
tistical mechanics for the study of criticality. Also one can hope to calculate the relevant 
quantities on the lattice. This goal is fundamental since lattice calculations remain prac- 
tically the only non-perturbative method in QFT that can be derived strictly from first 
principles. Thus, it is certainly worthwhile to continue the exploration of possibilities to 
obtain high-energy or small-xej cross sections from Euclidean quantities. 
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Appendix 

A Considerations concerning the analytic continua- 
tion of amplitudes in the z/ plane 

Here we will first discuss the analytic continuation of Tr(i^, Q^), Eq. ( |2.12| ), from real 
values of v into the half plane Imz/ > 0, where we have 



Imz/ > Im^z/2 + g2 > 0. (^1) 

Proof: the second inequality is obvious. The first inequality holds for imaginary v. By 
continuity, the inequality can only be violated if there exists some v in the upper half 
plane such that 



Im z/ = Im ^1^2 ^ Q2 _ (^2) 

Thus, for some a & M, 

g2 = j, + a. (0.3) 



2 



Squaring this equation, one finds 

Q2 = 2iya + , (0.4) 

which can not be fulfilled for non-zero Q and Imu > 0. Therefore, the first inequality in 
Eq. (0.1) holds everywhere in the upper half plane and our proof is complete. 

With Eq. (0.1) we find that the factor exp(ix°z/ — ix^^/u"^ + Q"^) in Eq. ( ^.121 ) de- 



creases exponentially for Imz/ > as x — * oo with |x |/x < 1. Thus the integral in 
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Eq. ( p.l2 ) should be well convergent and define an analytic function in v for Imz/ > 0. 



There is still the possibility of T^iv^Q^^ having a cut on the imaginary v axis for 
< Im z/ < Q due to the explicit factor ^/l7^~+'Q^ in Eq. ( |2.12| ) . But using Eq. it is 



easy to see that the values of Tr(z/, Q"^) when approaching the imaginary axis from both 
sides are equal. Then the "edge of the wedge" theorem (cf., e.g., [^) guarantees that 
Tr{i^, Q^) is also analytic for < Imi/ < Q. The point v = iQ could still be an isolated 
singularity but this is incompatible with the square root in Eq. ( ^^.121 ). This concludes 
our discussion of the analyticity of Tr(z/, Q'^) for Imi/ > 0. 



The discussion of the analyticity of T+(z/, Q^, /i) defined as in Eq. ( |2.14D but with irj 



replaced by z/ is analogous. The proof of Eq. (|2.15|) for < ?7 < Q is straightforward 
using Eq. (|2Tl| ). 

In the remainder of this appendix we will discuss the connection between the limits 
z/ — s> oo either on the real or the positive imaginary axis. We will outline one possibility 
of specifying the requirements for guaranteeing that the analytic continuation of the 
asymptotic behaviour of the amplitude from the real to the imaginary axis is possible. 

The retarded amplitude Tr{v,Q'^) is an analytic function of v with cuts along the 
positive and negative real axis (Fig. |[). Let us assume that there exists a function /(z/) 
analytic for Imz/ > such that f{v) ^ for Imz/ > 0, and 

is a bounded function, |-R(z/)| < const, for Imz/ > 0. (Since is kept fixed we suppress 
this argument in -R(z/) and /(z/).) Let us furthermore assume that 

\i^^R{r,) = C^. (06) 

Then the Phragmen-Lindelof theorem (see theorem 5.64 of ||2^) states that C+ = C_ 
and that -R(z/) approaches the same limit C+ along any ray v = r]exp{i(f)), < < vr, 
= const., 7] — > +00. Thus we have also on the imaginary axis 



lim R{i7]) = C+ . (0.7) 



Suppose now that from the study of Tr{iri, Q^), e.g., by lattice calculations, we can deduce 
the asymptotic behaviour for r] oo, construct a suitable function f{iri) and deduce the 
value of the constant C+ [Eq. (0.7)] with C+ ^ 0. 

Typically one would try to fit lattice data to functions reflecting our general expec- 
tations concerning the high-energy behaviour of amplitudes in QFT like /(z/) = (z/ — z/q)" 
or /(z/) = (z/ — z/o)"ln(z/ — z/q). Here a should be taken as a real constant and z/q as a 
constant with Imz/o < 0. With the assumptions specified in Eqs. (0.5) and (0.6) we get 
for the physical amplitude on the real axis: 

T,(z/,g2) — > C+f{u) for z/^oo. (08) 

To summarize: if we assume that the limit of T^(z/, Q^) for z/ —>■ oo, < argz/ < vr is 
governed by a suitable analytic function /(z/), then this function can be determined trom 
the study of Tr(z/, Q^) on the positive imaginary axis. 



21 



We have only outlined here the simplest assumptions making the analytic continua- 
tion of the high-energy limits on the real and imaginary axes possible. No doubt, using 
the methods of one can relax these assumptions and still get useful relations between 
the two limits. 



B The functional integrals for A4e cind A4m 

In this appendix we give the details of the derivation of Eq. ( [4.131 ). We start with 
Eq. (p.l3|) and using Eq. ( |4.3| ) rewrite it as 

ME(T,r) = ^ lim lim e^^f-^^^^UOl e'^^f-^f^^^f^'^''^ A(0) e'^^f-^^^^ff^'^ (B.l) 



X J(0) e-(^-°)^=//('^) J(0) e-(°-^»)^^//W A\0) e-(-»~<)^'=//W |0) ^ e-^^'r<)''^ . 

Here we have introduced further times r/, t'j with t[ < Ti < Tf < t'j- and the limits in 
Eq. (0.1) are to be taken in the order indicated there. 

Let 10) be eigenstates of the field operator $(a;) at x° = 0: 

«l>(f,O)|0) = 0(x)|0), (02) 

where 0(5;) are classical functions. We have 

I D<P\<P){<P\ = 1. (03) 



The basic relation of the path integral formalism discussed in many textbooks (see 



^0(2) I ^-ArH.ffir) |0{1)^ ^ J ^^(2) J ^3477(2) (x) (<i!.{2) (x)) +iAr (tt^ (x),<?i{2) (x))] 

+0{At') , (04) 

where the measure includes a factor l/(27r) for each 7r'^^^(x) integration. Here Tr^'^\x) is 
the classical momentum field and Tir the classical Hamilton density corresponding to 
Heffir): 

Hrin, 0) = ^nixf + (1 - r)7rix)dsm + ^ (V0(f))' + ^m'<P{xf + ^0(f)4 . (^5) 

Next we choose a grid on the time interval [r/, rj-]: r/ = r^^^ < r^^^ < r^^^ ■ ■ ■ < r^^^ = t'j, 
where also r^, 0, r, Tf should occur as intermediate points t^^\ We discuss the factor 
exp[— (tj — T-)Eo] which we rewrite using Eqs. (0.3) and (0.4) as: 

^-ir'^-r^Eo ^ (0|e~("/-"^)^="'('')|0) (B.6) 

= (Ole-^"''^'-"*^"'^)^^//^^) • • • e-("*''-"'°')^-//(") |0) 



AT 



n (O|0(^))(0(^)| e-(-'"'--^"-^^)^=//M |0(^-^)) 

j=0 
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N N 



g-(r;,-rOi?o _ J fl D(l)^^^ II D7i^^\0\(t)^^^) (B.7) 

At^O j=0 j = l 

N 

X Y\_e^ /d^xO)[7r{#J"',T(j))((/)(x(j'),T(j))-</)(^(j),TO-l)))+i(TO)-TW^ (^'^'^^ | O) 

i=i 

Here we set (p^^^x) = r^-^^) and ir^^^x) = 7i(x,t^^^). 

Now the integration over the vr fields in Eq. (^.7) can be performed since the integral 
is of Gaussian type. In the limit N oo, Ar ^ we get 

where CE,eff is given in Eq. ( [4.15| ) and C is a constant (which is infinite). Performing 
the same steps for the matrix element in the numerator in Eq. (0.1) and taking the limit 
t[ — > — oo, Tj — > +00 we get the expression for Jvi e{t, r) of Eq. ( [4.13| ). If we perform the 
analogous steps for the Minkowskian matrix element, Eq. ( |3.19| ), we get: 

MMit,r) = i lim e^(*/-*0A/2-i / Dcj) aitf)jit)jiO)a\u) I '''^''^'■^ff , (^9) 

Zj '■i^ J 

t j7 — *+00 

Z = J D0e*/^'^^*^>=-f-f , (glO) 



with 

and 



^M,eff = l{do<P{x) - (1 - r) dsmf - l{V<P{x)y - ^0(x)2 - A<^(x)^ (ill) 



a{tf) = / d^x (ido(j){x) - i{l - r) d^(t){x) + M0(x)) , (B.12) 
a\ti) = ( rf^f (-i9o0(x) +z(l -r)930(x) +M0(x)) . (§13) 
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